Reflectionless tunneling in ballistic normal-metal— superconductor junctions 
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We investigate the phenomenon of reflectionless tunneling in ballistic normal-metal- 
superconductor (NS) structures using a semiclassical formalism. It is shown that applied magnetic 
field and superconducting phase difference both impair the constructive interference leading to this 
effect, but in a qualitatively different way. This is manifested both in the conductance and in the shot 
noise properties of the system considered. Unlike diffusive systems, the features of the conductance 
are sharp and enable fine spatial control of the current, as well as single-channel manipulations. We 
discuss the possibility of conducting experiments in ballistic semiconductor-superconductor struc- 
tures with smooth interfaces and some of the phenomena specific to such structures that could 
be measured. A general criterion for the barrier at NS interfaces, though large, to be effectively 
transparent to pair current is obtained. 
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I. INTRODUCTION 

One of the most interesting phenomena in hybrid dif- 
fusive normal-metal-superconductor structures is reflec- 
tionless tunneling. This phenomenon manifests itself as 
a zero bias peak in the differential conductance of a dif- 
fusive normal metal slab connected to a superconductor 
via a tunnel barrier with low transmission probability T 
van Wees et al. ||] used a path integral picture to 
suggest and explain the effect of reflectionless tunneling. 
They show that the enhanced conductance at zero-bias 
is due to electron-hole coherence in trajectories that due 
to the disorder in the normal metal hit the barrier at 
the normal-mctal-superconductor (NS) interface many 
times. This results in the barrier being effectively trans- 
parent to pair current. 

In this paper we show that the phenomenon of reflec- 
tionless tunneling exists also in ballistic systems, the re- 
quirement being the existence of multiple reflections from 
the NS interface due to the geometry of the structure. As 
in diffusive systems, we find an enhanced NS conductance 
for zero-bias and zero magnetic field. We show that the 
magnetic field (iJ), finite energy and voltage, and super- 
conducting phase difference (<I> S ) impair the constructive 
interference leading to the enhanced NS conductance, 
but applying the superconducting phase difference has 
qualitatively different consequences than applying a fi- 
nite magnetic field or voltage. 

We show that the ballistic nature of the system gives 
rise to pronounced and delicate features, which are not 
averaged over as in the case of diffusive systems. This 
results in new measurable phenomena, such as sharp 
peaks in the NS conductance as new channels open, and 
quasiperiodicity of the conductance as a function of mag- 
netic field. We also demonstrate the possibility, specific 
to ballistic systems, to conduct detailed manipulations 
such as extracting out a single channel from a normal 
metal (semiconductor) waveguide or extracting the cur- 
rent at a given position along the waveguide. 

The ballistic regime in semiconductor-superconductor 



hybrid structures was investigated recently experi- 
mentally PHlC]. Unlike the case in normal- metal- 
superconductor structures, where sharp boundaries are 
made that enable specular reflection at the NS inter- 
face [II 12], in semiconductor-superconductor interfaces 



specular reflection is sacrificed for the purpose of lower- 
ing the barrier at the interface, thus increasing the An- 
dreev reflection probability. We here raise the possibil- 
ity to conduct experiments in ballistic semiconductor- 
superconductor structures with a sharp interface and a 
long elastic mean free path. Though indeed the transmis- 
sion probability of the barrier would then be small, the 
electron- hole coherence over long trajectories results in a 
large Andreev reflection probability, as we show below. 
Thus, one can have strong proximity while preserving 
the ballistic nature of the system. Other systems which 
seem favorable for the realization of ballistic NS struc- 
tures with specular reflection at the interface are the re- 
cently investigated organic molecular crystals Ju|[lJ|] • In 
these systems the NS transition could be realized by ap- 
plying a space dependent gate voltage. 

The paper is arranged as follows: In Sec. || we intro- 
duce the formalism and the structure we consider, ob- 
tain the expressions for the three-terminal conductances 
in terms of Rhe{N), the Andreev reflection probability of 
a trajectory that hits the interface N times, and calcu- 



late this probability for zero magnetic field. In Sec. Ill 
we show that for a short slab the NS conductance has 
sharp peaks as channels open. In Sees. ^ and [v] we cal- 
culate Rh e (N) and the linear conductances as function 
of H (fyj) and $ s in a similar SNS structure (0). In 
Sec. VI wc calculate the shot noise in both structures, as 



a function of H and <!>». In Sec. VII we consider diffusive 



systems and demonstrate the connection between the ef- 
fect of reflectionless tunneling in diffusive and ballistic 
systems. Throughout the paper we consider zero tem- 
perature and use the model in which the superconduct- 
ing order parameter A is constant in the superconductor 
and zero in the normal metal. 
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II. CONDUCTANCE OF A LONG NORMAL 
SLAB ATTACHED TO A SUPERCONDUCTOR 

A. Model 

We consider a ballistic normal-metal or semiconduc- 
tor slab between two normal reservoirs. The slab is sep- 
arated by an infinite barrier from a region denoted as 
vacuum, except in a region of length L, at which a su- 
perconductor is attached to the slab (Fig. [IJ) . At the NS 
interface the barrier is finite, with transmission probabil- 
ity r. The opening of the normal slab to the two normal 
reservoirs is taken to be adiabatic and the length of the 
slab between the reservoirs and the NS interface to be 
long enough such that channels are formed with homo- 
geneous distribution in the transverse direction . We 
also assume that the change from infinite barrier to fi- 
nite barrier of transmission T at the end points of the 
NS interface is not abrupt, but smeared over a length s 
such that Af <C s C W, where W is the width of the 
slab in the direction perpendicular to the interface. In 
this way the change is adiabatic but the smearing can 
be neglected in our calculations. We denote this struc- 
ture as a vacuum-normal-metal-superconductor (VNS) 
structure, as opposed to a similar structure with another 
superconductor attached symmetrically to the other side 
of the slab, which will be an denoted SNS structure. 

The superconductor is connected to a third reservoir 
except when explicitly mentioned otherwise. We con- 
sider the case where the electrochemical potentials of the 
right and superconducting reservoirs are equal, and the 
left reservoir is biased by an infinitesimal voltage, and 
calculate the three-terminal linear conductances of the 
system. Previous works concerning similar structures 
]Tg| |l9t considered the NS interface either as fully trans- 
parent or concentrated on effects of channel mixing due 
to the roughness of the barrier when it exists. We con- 
sider the NS interfaces to have a smooth barrier, so that 
normal reflection is specular and the Andreev reflected 
hole retraces the electron's trajectory. We assume spec- 
ular reflection from the VN interface as well. 

Our model is two dimensional. While assuming Xp <C 
W, therefore having many channels, we assume for sim- 
plicity that the thickness of the slab (the third dimension) 
is small, having one transverse mode in this direction. 
The generalization of our treatment to thicker slabs is 
trivial. 

We use a semiclassical formalism and consider the 
propagation of electrons in each channel to be described 
by their classical deterministi c trajectory PJ20|]. For each 

channel j we define kj\\ = \J 2mEp /Ti 2 — j 2 ir 2 /W 2 and 

calculate the angle 6j — ta,n~ 1 {jir/ (kj\\W)] between the 
classical trajectory of an electron in this channel and the 
NS interface. We consider an electron entering the nor- 
mal slab from the left reservoir, approaching the region 
of the slab with the NS interface ( "NS region" ) at a given 
distance from the NS interface and angle 9j with respect 



to it. If the electron is only normally reflected from the 
NS interface, it follows a certain trajectory in the slab 
until exiting it to the right reservoir after hitting the NS 
interface N times ("N trajectory"). Due to the finite 
Andreev reflection amplitude at each point it hits the 
NS interface, the electron has a probability Rhe(N) to 
be reflected as a hole to the left reservoir. In this model, 
due to the interfaces being parallel and smooth, there is 
zero probability for an electron to be reflected back to 
the left reservoir or to be transmitted as a hole to the 
right reservoir. Therefore, Rh e (N) +T ee (N) — 1, where 
T ee (N) is the probability of an electron coming from the 
left reservoir to be transmitted as an electron to the right 
reservoir. 




FIG. 1. Vacuum-ballistic normal- metal-superconductor 
junction with a barrier at the NS interface. Each time the 
particle hits the NS interface it can be reflected either nor- 
mally or in an Andreev process. 1 and 2 are normal reservoirs 
and 3 is a superconducting reservoir. 6 is the angle of inci- 
dence. Filled (empty) arrows designate electrons (holes). In 
Sec. |v| a similar structure is considered, with a second super- 
conductor attached to the other (lower) side of the slab, so 
the structure has up-down symmetry. 

For each open channel in the slab the number of times 
a trajectory hits the NS interface is either Nj or Nj + 1, 
where Nj equals the integer part of Lt&n9j/(2W), with 
L being the length of the NS interface. The fraction of 
trajectories in channel j that hit the NS interface Nj + 1 
times is given by Pj = L tan 9j / (2 W) — Nj . The Andreev 
reflection probability of an electron in channel j is then 
given by 

R j he = Pj R he (Nj + 1) + (1 - P j)R he (Nj) . (1) 

We define by I\ , Ii and I3 the currents emerging from 
the left terminal, right terminal, and superconducting 
terminal, respectively. Due to current conservation I\ — 
—I2 — ^3- We then define the NN, NS, and total linear 
conductances of the system as 

= H v = ¥ £ e(fc 'i )(1 - RL) ' (2) 

G 3i = - to^jEe^, (3) 
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and 



G T = lim | = G 21 + G 31 = ^ Yl e (*y||)( 1 + O > 



where 0(x) is the Heaviside theta function. 



(4) 



B. Andreev reflection probability of an N trajectory 

The calculation of the conductances is therefore re- 
duced to the calculation of Rh e {N) = |r/ le (7V)| 2 where 
r he{N) is the corresponding amplitude. For a single hit 
at the NS boundary, we denote by r^e (feh) the ampli- 
tude for an electron (hole) to be Andreev reflected and 
by r ee (rhh) the amplitude for an electron (hole) to be 
normally reflected. By dividing an N trajectory to an 
N — 1 trajectory and a 1 trajectory, we obtain a recur- 
sion formula 



rhe(N) = r he + r ee r he (N - l)r hh -\ 
reeThe{N - l)r eh r he (N - l)r hh + . 

r e erhe(N - l)r hh 



= r hl 



1 - r eh r he {N - 1) 



(5) 



Using the relations (which are exact at Ep) 



ee — ' hh' 



r e h = fhei and \r e h\ 2 + |r ee | 2 = 1, we assume, and then 
show by induction, that r} le {N) is imaginary for all N 
and can be written as 



r he {N) = i 



1^(^-1)1 + 1^1 



(6) 



l + |r he ||r fce (JV-l)| ■ 
The solution of this equation is given by 

r he (N) = itanh[JVtanh _1 (|r he |)] . (7) 
For a barrier with small transmission probability we find 
Rhe(N) w tanh 2 (iVr) « tanh 2 (AT/2) , (8) 

where we define r — |r^ e | = T/(2 — T). 

Using Eq. (^) to obtain the values of Rh e (N) for all 
the channel-dependent Nj and Nj + 1 in the conductance 
formulas [Eqs. (§)-(§], we obtain the linear conductances 
G21, G31, and Gt- In this paper we are interested in the 
case were r< 1 and, therefore, in some of the formulas, 
and in the qualitative discussions, we take this limit. 

Before considering further the conductances of the sys- 
tem, we would like to dwell on the physical aspects of 
Eq. (||). This formula reflects the essence of the physics 
behind "reflectionless tunneling." It states that electrons 
in trajectories that hit the NS interface N 3> 1/r times 
are Andreev reflected with probability close to unity, even 
though r< 1, thus making a barrier having a low trans- 
mission coefficient effectively transparent to pair current. 
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FIG. 2. Vacuum-ballistic normal-metal-superconductor 
junction where a part of the normal metal is removed (only 
the relevant region is shown). 

This is a result of electron-hole coherence in the nor- 
mal metal. For an incoming electron, the different paths 
resulting in a hole returning to the reservoir interfere 
constructively, while the different paths resulting in an 
electron transmitted to the right reservoir interfere de- 
structively. The constructive interference for a returned 
hole competes with the small Andreev amplitude at each 
encounter with the interface, which is proportional to 
T, and therefore R he {N) « 1 only for N > r -1 . This 
means that for channels in which tan# 3> 2W/{YL) the 
barrier at the NS interface is not effective. In fact, if 
one considers a system in which the superconductor is 
floating and the above condition is fulfilled for all the 
channels, one can show that the current between the left 
and right reservoirs flows inside the superconductor and 
a part from the middle of the normal slab can be taken 
out (see Fig. ||) without affecting the conductance of the 
system. 

Equations (0) and (@) are far more general than the 
above model and hold in any case where an electron in 
the normal metal can hit the NS interface more than once 
before electron-hole coherence is lost. This is true in var- 
ious geometries in ballistic systems, and also in diffusive 
systems, which are considered in Sec. VTJ. In all these 
cases, Eq. (g) results in a criterion for the effectiveness of 
a barrier with small transmission probability: Consider a 
physical property which is determined by a certain set of 
trajectories; the criterion for the barrier at the NS inter- 
face not to be effective is that most of these trajectories 
hit the interface more than T" 1 times before electron-hole 
coherence is lost. In Sec. VII we will show how this gen- 
eral criterion reduces, in diffusive systems, to the known 
conditions for the barrier, though high (r <C 1), not to 
affect the conductance and the density of states of a dif- 
fusive NIS junction. 



C. Comparison to an incoherent structure 

In order to show that Eq. (||) is a result of constructive 
interference, which is due to electron-hole coherence, we 
compare our result to the case where, due to strong de- 



3 



phasing, there is no electron-hole coherence, i.e., where 
the phase between two consecutive hits of the interface 
is lost. In this case the problem is reduced to a ran- 
dom walk problem, with forward-backward asymmetry. 
At each hit at the NS interface the electron (hole) has a 
probability T 2 <C 1 to be Andreev reflected, in which case 
the direction of propagation is reversed and the proba- 
bility to move forward is 1 — T 2 . The size of the step is 
channel dependent and is given by dj = 2W cot 6j, the 
distance between two consecutive points a trajectory in 
channel j hits the interface. This gives a mean free path 
of lj ~ dj/T 2 . For L <C lj the Andreev reflection prob- 
ability is Rhe = L/lj, and for L ^> lj it is 1 — lj/L (the 
probability for a transmitted electron is lj / L) . Therefore 
lj is the saturation length, beyond which the Andreev 
reflection probability is close to unity. On the contrary, 
in the case of coherent scattering, using Eq. (||) and the 
relation Nj — L/dj, we find that the saturation length is 
dj/T = IjT. Due to the scattering being coherent, it is 
smaller by T compared to the noncoherent case. More- 
over, for short slabs, L <C IjT, R he ss L 2 /(Aljdj), larger 
by L/(4dj) than the noncoherent case. For long slabs 
(L ^> IjT ) one obtains Rhe ~ 1 — exp(— Lj ^/2ljdj), 
and the probability for a transmitted electron is expo- 
nentially small, and not linear in lj/L as in the nonco- 
herent case. The difference between the two cases is most 
notable for slabs with intermediate lengths between the 
two saturation lengths, IjT <C L <C lj, which corresponds 
to 1/r < N < 1/r 2 . Without coherence R he < 1, and 
with coherence Rhe ~ 1- 

It is instructive to compare Eq. (||) to a similar sys- 
tem, in which the superconductor is not attached to the 
slab on its side, but part of the slab itself, of length L, 
is superconducting. In this case and assuming no barrier 
at the NS interfaces, the Andreev reflection probability 
of an incoming electron is given by tauh 2 (AL/2hvF) = 
tanh 2 (L/2£ s ) pi] . Here it is the ballistic supercon- 
ducting coherence length, which is the length scale for 
pairing. In our system, Eq. (^J) can be written as 

R he {L) «tanh 2 [ J Lr/(2^-)] , (9) 

with the saturation length dj/T as the length scale for 
pairing in the normal slab due to the proximity to the 
attached superconductor. 

III. CHANNEL OPENING 

Using Eqs. @ and ^ the NS linear conductance can 
be calculated as function of Ep, T, L, and W . We now 
concentrate on a special case of these parameters, which 
results in sharp resonances of the NS linear conductance 
as function of the Fermi energy. While in all the other 
cases considered in this paper we are interested in the 
emergence of constructive interference that leads to en- 
hanced Andreev reflection and therefore consider cases 
in which electrons in at least some of the channels hit 



the interface N ^ 1/T times, we are now interested in 
a different limit, in which TL^Jkp /W <C 1. The generic 
behavior in such a structure would be that the current 
flow to the superconductor is small, since the number 
of times an electron in any transverse channel hits the 
barrier is smaller than 1/r. However, if we change the 
Fermi energy (e.g. by a back gate) such that the channel 
of the highest transverse mode has just opened, then the 
trajectory of a particle in this channel is almost perpen- 
dicular to the interface. As a result, the number of times 
a particle in this channel hits the interface is much larger 
than 1 /r and the contribution of this channel to the NS 
conductance is significant and given by 
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FIG. 3. (a) The conductance (in units of 4e 2 /h) between 
the left normal reservoir and the superconductor is plotted as 
function of E F [in units of h 2 /{mW 2 )] for LT/(4:W) = 0.1. 
(b) Enlarging the first peak we see that the conductance at 
the peak is unity, and the width of the peak is approximately 
[Lr/(4M/)] 2 = 0.01 of the value of E F at the peak. 
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In this equation we assumed that r < 1 and neglected 
the fact that a fraction of the trajectories hit the barrier 
Nj+1 and not Nj times, since as the channel opens Nj ^ 
1. Defining ej = h 2 (k jl{ ) 2 /(2m) and e = [LY/(AW)] 2 E F 
we find that for < ej <C e 



G sing = 4e_ tanh 2 
61 h 




4e" 



(11) 



While in a normal quantum point contact connected to 
a superconductor in series the linear NS conductance as 
function of the Fermi energy in the slab would show steps 
p2j, in our case, where a superconductor is attached to 
the point contact on its side, with a barrier at the inter- 
face, the NS conductance as a function of Ep has sharp 
peaks at the energies where channels open. The magni- 
tude of these peaks is 4 times the quantum conductance, 
and the scale of their ene rgy wid th is e. With the con- 
dition given above, TLy/kp /W <C 1, these peaks are 
narrower than the energy difference between the opening 
of adjacent channels. For L m W this condition reduces 
to T < y/Xp/W. Under the semiclassical approxima- 
tion we make the conductance peaks are nonanalytical as 
a function of Fermi energy at energies where transverse 
channels are opened, as can be seen both in Eq. ([ll]) and 
Fig. [| This is due to the fact that the number of times 
an electron hits the interface diverges as Cj — > 0. These 
nonanalyticities are a consequence of the semiclassical 
model, and one has to take into account that the validity 
of the semiclassical approximation is limited by the con- 
dition kj\\ 1/L*, where L* is the range of the potential 
variation. We estimate L* by fc^/Vfc^ = W * s / '(Xp\/T) , 
taking into account the variation of k±_ with the spatial 
variation of T. The condition kj 1/L* is equivalent 
to €j » [TXp/(W 2 s 2 )]Ep, which is consistent with the 
condition e~ <C e given that T ^> X F /(L 2 s 2 ). One can 
therefore expect that with this condition fulfilled, going 
beyond the semiclassical approximation would smoothen 
the above nonanalyticities, but will not alter the other 
features of the peaks (height and width) . For €j 3> e the 
contribution of the jth channel to the NS conductance is 
proportional to 1 /tj and is given by 



^31 - "IT — 



(12) 



In this section we described the effect of channel open- 
ing on the NS conductance of the system. The behavior 
of Gp is similar, only the peaks at the energies where 
channels open are half the magnitude and are on top of 
the step function of magnitude 2e 2 /h (since Gp is simi- 
lar to G 3 i, only 2R j he is replaced by 1 + B? he ), The NN 
conductance is given by the complementary of half the 
NS conductance to a step function {2R 3 he is replaced by 
1 — R° he )- In the next sections we mostly consider G31, 
and analogies to Gp and G21 can be made in a similar 
way. 



IV. MAGNETOCONDUCTANCE OF A LONG 
NORMAL SLAB ATTACHED TO A 
SUPERCONDUCTOR 

In this section we consider the same VNS structure 



as in Sees. |II| and III with a magnetic field applied per- 
pendicular to the slab. We investigate the effect of the 
magnetic field on the transmission probability of an ./V 
trajectory, as well as on the linear conductances in the 
system. Magnetic field penetration into the supercon- 
ductor is neglected. 

We consider H <C <&q/(XfW) where $0 is the flux 
quantum. Under this condition the curving of the trajec- 
tories of the particles in the normal slab can be neglected 
p3| . The number of times, Nj (or Nj + 1), a trajectory 
of an electron in channel j hits the NS interface stays un- 
changed, and so do Eqs. (||)-(||), except that the Andreev 
reflection probabilities Rhe(N) now depend also on 
the phase acquired by an electron and a hole moving in 
opposite directions between two consecutive points the 
trajectory hits the NS interface ("trajectory section"). 
This phase is given by $h — 4:irHA/&o, where A is the 
area of the triangle enclosed by the trajectory section and 
the interface. 

We now turn to the calculation of Rh e (N,$H)- Re- 
peating the same procedure leading to Eq. (|b|), but keep- 
ing track of the phase introduced by the magnetic field, 
we obtain 

r ee r he (N - l,$ H )e i ' s>H r hh , . 
l-r eh r he {N- l,$ H )e l *H 

In order to obtain an explicit formula for r^ e (iV, $#) it 
is useful to define 



In 



(14) 



where 0n and 7jv are also <&h dependent. Inserting this 
definition for N and N — 1 into Eq. ( |l3| ) we obtain the 
matrix equation 



(3n 
In 







1 


= 




1 




f e i*H 


1 


= 




1 



JV-1 



(15) 



Here we used the fact that at the Fermi energy r^ e is 
imaginary. /3/v and 7^ are in principal defined up to (the 
same) multiplication constant, which we dictate by the 
choice P\ = 71 = I . Diagonalizing the matrix and taking 
it to the power N — 1 we obtain 



r he (N, <f> H ) 



(16) 



-i* H /2 



i sin ($ff/2) + v^coth 



N tanlT 



cos($tf/2) 
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Rh 



where b = r 2 — sin 2 (<f>#/2). Using the fact that the sec- 
ond term in the denominator is always real (also when b 
is negative) we obtain 



R he (N, $ H ) 



(17) 



sin 2 ($ H /2) + &coth 2 



N tanh 



Vb 



cos($ ff /2) 



which, for b < 0, can be written as 
R he (N,$ H ) = 



(18) 



sin 2 ($ij/2) + (-6) cot 2 



iV tan- 



cos ($ff/2) 



These formulas hold for any JV, 
now consider the case of r <C 1 
sin 2 (<E>#/2) > r 2 (which holds 
when r <C 1), we see from Eq. (16 



r(r), and ff. We 
and iVr > 1. For 
for most values of 
_J that ifo e (iV, $h) < 1. 
In the opposite limit, of sin 2 (<J>#/2) <C r 2 , one obtains 
from Eq. Q that R he (N,$ H ) w 1. This leads to the 
conclusion that the Andreev reflection probability from 
an iV trajectory is small for almost all values of perpen- 
dicular magnetic field, except those special values that 
result in \$ H - 2kn\ < r (Fig. |). Between every two 
such peaks the function oscillates, having N — 2 smaller 
peaks (and N — 1 nodes). These peaks (nodes) corre- 
spond to integer (half integer) flux quanta through an 
area of an integer number of triangles in the trajectory 
of the specific channel. 

The magnetic field not only impairs the constructive 
interference leading to large Andreev reflection at zero 
field, but causes destructive interference. This can be 
seen by considering <& H = 3>o/2. Then, R he (N) = 
for even TV and for odd N it equals r 2 , the Andreev 
reflection probability from a single hit. Therefore, for 
any given cha nnel Ri e is of order r 2 . In the cases dis- 
cussed in Sec. II C, of no interference and of constructive 



interference, there were (different) saturation lengths be- 
yond which the Andreev reflection probability was close 
to unity. Here, however, due to destructive interference, 
there is no such length scale and Andreev reflection is 
small ( r 2 ) for any length of NS interface. This is true 
also when the destructive interference is due to a super- 
conducting phase difference in an SNS structure, as is 
discussed in Sec. jvj. 

Inserting Eq. (|l7|) in the conductance formulas 
[Eqs. (d)-(||)] we obtain the linear conductances as 
function of magnetic field for any junction parameters 
(r, L, W, Ep). If the parameters of the structure are such 
that TL/W ^S> 1, then at zero magnetic field the NS di- 
mensionless conductance is much larger than unity, since 
there are many channels for which NT ^s> 1. We choose 
such a case and plot in Fig. || the NS dimensionless con- 
ductance as a function of H . At H = the conductance 
has a sharp peak, of mag- 
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FIG. 4. R he (N, $ H ) as obtained from Eq. (0) for (a) 
N = 10 trajectory and (b) JV = 50 trajectory in a VNS sys- 
tem (Fig. M) for a barrier transmission probability F = 0.2. 
Notice the narrow large peaks periodic in and the small 
oscillations between each such peaks, having in (a) TV — 1 = 9 
nodes. In (b) the magnitude at the high peaks is approxi- 
mately unity, and the oscillations between them are hardly 
visible. 

nitude of the order of the number of channels and width 
of order r. As H is increased, the constructive interfer- 
ence leading to the enhanced Andreev reflection is de- 
stroyed in one channel after the other and the NS con- 
ductance becomes small. However, the conductance of 
each channel is periodic in H, with a period given by 
the area of the triangle between a trajectory section in 
this channel and the interface. This quasiperiodicity is 
reflected in the peak spectrum of the NS conductance 
as function of H (Fig. ||). Periods of larger H reflect 
channels with a smaller triangle, which corresponds to 
a trajectory of larger N, and therefore [see Eq. (||)] the 
peak heights are larger. Using this quasiperiodicity one 
can obtain "magnetic switching." By choosing the mag- 
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FIG. 5. The conductance (in units of 4e 2 //i) between the 
left normal reservoir and the superconductor as function of 
H (in units of <$> /W 2 ). The plot is given for L/W = 50, 
T = 0.1 and 2mE F W 2 /Ti 2 = 1000, which corresponds to 10 
open channels. Peaks higher than unity are a result of over- 
lapping resonances of 2 or more channels. Apparent periods 
of H are: 0.29, 0.37, 0.48, 0.65. At H = the peak is signifi- 
cantly higher than the others (G31 = 5, not shown). 

netic field such that \& H = 2fc7r| for one channel only, 
one can remove only electrons propagating in this chan- 
nel from the normal slab to the superconductor and have 
the electrons in all the other channels propagate to the 
right reservoir with probability close to unity. 

Though the results in this section were given for the 
linear conductance at finite magnetic field, it is straight- 
forward to generalize our results to be valid for finite 
subgap voltage, thus obtaining the differential conduc- 
tance as a function of voltage. One can also incorporate 
a constant phase gradient V0 in the superconductor in 



parallel to the NS interface, generated by a constant su- 
percurrent. The Andreev reflection amplitude from an 
N trajectory would then be given by Eq. ([IB]) with <E># 
replaced by 

$ tP = $jr + V0d,-+#f(e). (19) 

Here 3>J r (e) = (k~^» — k~n)W/ sin 9 j+axg[r ee rhh] is the rela- 
tive electron-hole phase accumulated due to finite energy 
in one triangle, where fc^jj = kj\\(Ep — > Ep ± e). Note 
that there is a complete analogy between applying a per- 
pendicular magnetic field and constant gradient of the su- 
perconducting phase, with the relation V$ = 2ttHW/^q. 
For H = and V$ = 0, we obtain a zero-bias peak in 
the differential NS conductance as a function of voltage, 
similar to the low-iJ behavior of the NS conductance as 
shown in Fig. |5|. 

V. CONDUCTANCE PARALLEL TO THE 
INTERFACE IN AN SNS SYSTEM 

We now consider a system in which a second supercon- 
ductor is attached symmetrically to the other side of the 
slab, so the structure has up-down symmetry. We con- 
sider the case in which the two barriers between the nor- 
mal slab and the superconductors have the same trans- 
mission probability T. Then, one can apply the same ap- 
proach we used in the previous sections, only count the 
number of hits of each trajectory at both interfaces. The 
calculation of Rh e (N, $ s ) is given in Appendix |A|. It is 
done in the same spirit as the calculation of Rh e (N, <&#), 
but is more elaborate since one has to distinguish be- 
tween even and odd times a trajectory hits the interface, 
and the recursion relations are more complicated. The 
result for even TV is 



Rhe(2N,$ s 



z+(l-z) coth 2 



N tanh" 



2rcos($ s /2V(l-z) 
1 + r 2 cos $, 



(20) 



where z — r 2 sin 2 ($ s /2). The result for odd N is similar 
and is given in Appendix |A|. These results are even in 
$ s and, therefore, the same for a trajectory hitting first 
either of the two superconductors. For r <C 1 we ob- 
tain R he (2N,$ s ) = tanh 2 [2Nr cos($ s /2)]. For Nr > 1 
we see that R he (2N,$ s ) « 1 unless |$ s - (2k + < 
w/(Nr), while R he (2N, $ 8 ) = for $ s = (2k + 1)tt. As 
a result, in this limit, the Andreev reflection probability 
from a 2N trajectory as function of <1> S has sharp dips 
near $ s = (2k + 1)tt of an approximate width of l/(2Nr) 
(Fig. ^|) . The transmission probability of electrons to the 
right reservoir is given by 1 — Rh e (2N, <£> s ) and has sharp 
resonant peaks, which indicates that there is a transverse 



Andreev level shifted to Ep at <i> s = (2k + 1)tt, similar to 
the case in standard SNS junctions |24| . As in the case 
of perpendicular magnetic field, Rh e (2N, Q s ) has a max- 
imum at $s = (or multiples of 27r), but there are two 
major differences between the dependence of Rhe(N) on 
H and on <& s : (i) In a period of 27r Rh e (N, $ s ) has one 
minima, where Rh e (N,<&H) has N minima, (ii) while 
Rhe(N^n) exhibits sharp peaks near $h — 2irn, and 
for most values of magnetic field (generically) construc- 
tive interference is lost, and Andreev reflection is small, 
the situation for Rh e (N, <fr s ) is opposite. It is close to 
unity for most values of phase difference, and exhibits 
sharp dips near <f> s = ir + 2wk. 
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FIG. 6. R he {N,$ s ) in an SNS structure for N = 50 and 
F = 0.2. Note the narrow dips, in comparison to the narrow 
peaks in Fig. [j|(b), which is drawn for the same N and T, as 
a function of H . 

These differences can be understood by examining the 
two mechanisms of the destruction of the constructive in- 
terference between different paths of the same trajectory 
that result in a hole returning to the left reservoir. In- 
deed, the phase difference between a hole resulting from 
an Andreev reflection at the first hit of the NS interface 
and a hole resulting from an Andreev reflection at the 
second hit of the NS interface is similar in both cases 
(&h and <& s ), but the phase difference between this hole 
and a hole resulting from an Andreev reflection at the 
iVth hit of the NS interface is very different for the two 
cases. It is (TV— 1)&h for the case of a magnetic field and 
Q s or (depending if ./V is even or odd) for the case of su- 
perconducting phase difference. This introduces a large 
amplification factor in the electron-hole phase difference 
introduced by the magnetic field compared to that in- 
troduced by the superconducting phase difference. As a 
result, the magnetic field is far more efficient in destroy- 
ing the constructive interference leading to the enhanced 
Andreev reflection. 

The linear conductances as function of <E> S are cal- 
culated by inserting the results fo r th e Andreev reflec- 
tion probabilities [Eqs. ( |20| ) and ( A1C| )] in the conduc- 



Nj in the VNS structure 



tance formulas [Eqs. 
is now replaced by Nj which equals the integer part of 
Lta,n9j [W and pj is replaced by pj — LtaxiOj/W — Nj 



In Fig. [7| we plot (solid line) the NS conductance as a 
function of the superconducting phase difference for a 
system with the same parameters as the one in Fig. [| for 
comparison. The NS conductance for a similar structure 
with a larger barrier transmission probability (r = 0.25) 
is also plotted (dashed line) to demonstrate the narrow- 
ing of the width of the dips as AT grows. Here we see 
another marked difference between applying a perpen- 
dicular magnetic field in the VNS structure and a phase 
difference in the SNS structure. In the case of the applied 
magnetic field there is a large peak at H — 0, to which 




FIG. 7. G31 (in units of 4e 2 /h) as function of $ s . The solid 
curve is given for the same parameters as in Fig. [B| to enable 
comparison. In contrast with the case of applied magnetic 
field, the conductance is periodic in <& s , has narrow dips at 
odd multiples of n, and the oscillations as function of 3> s are 
of the order of the full conductance (giant). The dashed line 
is plotted for the same parameters, except here F = 0.25. The 
larger conductance and the narrower dips are both a result of 
NT being larger for each channel. 

all the channels contribute due to the constructive inter- 
ference, but the periodic peaks at higher fields appear for 
each channel at a different H . If the parameters of the 
junction are such that the separations between conduc- 
tance peaks are smaller than their width, the result will 
be a smooth oscillatory behavior of the conductance as 
function of H. On the other hand, in the case of SNS 
structure, the dips at all <& s = (2k + 1)tt are common to 
all the channels, and therefore the conductance oscilla- 
tions as function of $ s show sharp features of magnitude 
of the order of the total conductance. 

Recently, Petrashov et al. measured large conductance 
oscillations as a function of the magnetic field and super- 
conducting phase difference in a normal slab connected 
to "superconducting mirrors" p^ , |2^ |. Our results can- 
not be directly applied to the experimental structures 
studied by Petrashov et al., since in the experiment the 
structures are different, the superconductor is floating, 
and there is finite scattering in the normal slab. How- 
ever, some features appear to be more general and exist 
both in the experimental results and in our calculations. 
These are the much larger magnitude and sharpness of 
the oscillations as a function of superconducting phase 
difference compared to the oscillations as a function of 
magnetic field. 

To conclude this section we now apply the results ob- 
tained above to show how one can get controlled current 
withdrawal from an electronic waveguide. Here we use 
the fact that as long as the phase difference between the 
superconductors is it, the electrons move in the slab as 
in a waveguide. By replacing the bottom superconduc- 
tor with a series of superconductors (Fig. IsJ) , each with a 
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controllable phase <pi and interface length with the slab, 
L, we can create a "switch" in which we can control the 
location where the current is drawn. We set ef>^ui = tt 
and 4% = 0. For all i ^ i an incoming electron from the 
left in an angle with tan 9 ^> (2W/LT) will be normally 
reflected as in a waveguide. However, when it reaches the 
ith superconductor Andreev reflection occurs, adding a 
Cooper pair to the superconductor. One can therefore 
inject a current from the left reservoir and draw it at any 
one of the superconducting slabs. 
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FIG. 8. SNS junction where the bottom superconductor is 
divided to pieces with controllable phase of the order param- 
eter (only the relevant region is shown). 



VI. SHOT NOISE 

In this section we calculate the shot noise as function 
of H in the VNS structure and $ s in the SNS structure, 
and show that the differences between applying magnetic 
field and superconducting phase difference are reflected 
remarkably in the shot noise properties of the systems. 

We define the quantities P w = 2 f™^ dt(A/ i (i)AJ i /(Q)) 
(where are normal terminal indices), which give both 
the shot noise and the cross correlators between current 
fluctuations at the two normal terminals. Anantram and 



Datta 27 considered the case where an arbitrary number 
of superconducting and normal terminals exist, with the 
restriction that the chemical potential in all the super- 
conducting terminals is the same, and obtained general 
equations for the current correlators. Using their equa- 
tions for our system we obtain 



Pu — Pt. 



Pv 



P( 



T}3 



), (21) 



where Pq — 2e| V\ (2e 2 /h). This formula is applicable 
for both the VNS and SNS structures, and the specific 
parameters of the junction as well as the H and <E> S de- 
pendence enter only into R 3 he . Notice the full positive 
noise correlations between the two normal terminals [^8) , 
which is a result of zero normal reflection to the same 
reservoir and Andreev transmission to the other reser- 
voir in our model. 

Due to the dependence of the shot noise on the func- 
tions R 3 he (l — Rl e )) it shows peaks at points where the 
Andreev reflection amplitude is neither close to zero or 
unity. For the SNS structure we consider, this results in 




-37T-27T-7T 7T 2 71 3 71 

FIG. 9. Pu in units of Pq as function of $ s for the geome- 
try where the normal slab is attached to two superconductors. 
The parameters of the system are the same as those in Fig. ^ 
(dashed line). 

a sharp feature near the values of $ s corresponding to 
dips in the NS conductance [<E> S = (2k + l)ir], which are 
common to all the channels. The form of the sharp fea- 
ture is two double peaks separated by a very sharp dip 
p9[ , as can be seen in Fig. ^. As function of H, in the 
case of one channel, a sharp feature appears near each 
value of H corresponding to a peak in the conductance. 
However, these points are channel dependent, and as was 
the case for the conductance, the presence of many chan- 
nels smears the sharp features as function of H . 

Our results for the noise are easily generalized to finite 
subgap energy (differential shot noise as function of bias 
voltage) and gradient of the superconductor phase in the 
same manner discussed at the end of Sec. 



VII. DIFFUSIVE NS JUNCTIONS 



The semiclassical approach, introduced by van Wees et 
al. to explain the phenomenon of reflectionlcss tunneling 
in diffusive NS junctions || was used to analyze numeri- 
cally other experimental results as well (see, e.g., Refs. || 
and H ) . This approach has proved useful in obtaining a 
qualitative understanding of the physical phenomena in 
various geometries which make the use of the standard 
methods (quasiclassical formalism, BdG equations) diffi- 
cult. As Eq. (||) applies to any NS system, ballistic or 
diffusive, it can be used to obtain analytical results within 
this approach. In this section we apply the semiclassical 
formalism to treat both the phenomenon of reflectionless 
tunneling and the reduction of the local density of states 
(DOS) across a diffusive NS junction. We show that both 
of these phenomena are a result of the large transparency 
of the barrier to pair current (although r <C 1 and under 
the co nditi ons given by the general criterion at the end 
of Sec. IIB) and thus stem from the same physical effect. 








N' 



FIG. 10. Geometry of the model 
jectory with N=2. 



an example of a tra- 



We also demonstrate the connection between the effect 
of reflectionless tunneling in ballistic systems discussed 
in this paper, to the one in diffusive systems. 

Since we consider the particle's trajectory in the nor- 
mal metal to be deterministic, our approach can be ex- 
pected to give correct quantitative results [Eqs. (|2^) and 
(p3|)] when the scattering potential varies slowly on a 
scale of a wavelength j|,[30|, and the sample is short 
enough such that classical dynamics will not develop 
phase space structures on scales smaller than H |n) . For 
short range disorder, this approach is expected to give 
results which are qualitatively correct ||. 

Following the treatment of Ref. || , but using our ana- 
lytical result for the Andreev reflection probability from 
an N trajectory, Eq. (^), we obtain the linear conduc- 
tance of a normal slab connected via a barrier to a su- 
perconducting reservoir (see Fig. [jl]), which for T <C 1 is 
given by 



G(V -»■ 0, H = 0) = -j- J2 t2 (! - T ) N ^ tanh 2 (AT/2). 

(22) 



N=0 



Here n is the number of channels and T is the average 
transmission probability of the normal slab (the mean 
free path divided by its length). In the two limits where 
f « T and f » T Eq. (||) reduces to 



2e 2 n T 2 



G(V -*0,H = 0) 



(T « T) , 



(r»T) 



(23) 



Unlike ballistic systems, in this case multiple reflec- 
tions from the interface are enabled by the disorder. For 
small disorder r<T, the conductance is proportional to 
r 2 , reflecting the fact that Andreev reflection is a two- 
particle process. Already in this limit the disorder in- 
creases the conductance by a factor of 1/T. For T<T 
the disorder is large enough to generate, with high prob- 
ability, trajectories with N ^> 1/T, as we show below, 



and therefore the barrier is not effective and the conduc- 
tance is linear in T. The conductance has a maximum 
for T « T, where G « (2e 2 n/h)T w (2e 2 n/h)T. 

Equation (||) differs in the r > T, T -> limit by 
a factor of 2 from the analogous formula obtained by 
Beenakker et al. |32| for short-range disorder. A de- 
tailed discussion of how this discrepancy results from the 
different assumptions in the two models is given in Ap- 
pendix [b| 

As the essence of the effect of reflectionless tunneling 
is the fact that the barrier, though high, is transparent 
to pair current, the condition for the barrier to be inef- 
fective was discussed for this phenomenon |53],[54| as well 
as for other phenomena, as the reduction of the DOS on 
the normal side of an N-insulator-S (NIS) semi-infinite 
junction |33j]. We now show, usi ng ra ndom walk theory, 



that the criterion stated in Sec. II B , for the barrier to 
be ineffective, can be reduced in the diffusive case to the 
different known conditions for each phenomenon. 

In Appendix ^ it is shown that the typical length of a 
diffusive trajectory between N consecutive times it hits 
the barrier is Ln w N 2 l n , where l n is the elastic mean 
free path in the normal metal (interestingly, there is no 
average length for an N trajectory; see Appendix |^). 
This is also the order of magnitude of the length of the 
longest loop in such a trajectory (a loop is a part of a 
trajectory between two consecutive points it hits the in- 
terface). Using this result for Ln and since large contri- 
butions to Andreev reflection arise from trajectories that 
hit the interface N ^> times before losing electron- 
hole coherence (]|), only when coherent trajectories with 
total lengths larger than Lr = l n /F 2 occur with high 
probability will the barrier not be effective. This requires 
the electrons and holes to be coherent over a distance 
\/Ly * l n = 'n/r from the interface. Therefore, the gen- 
eral condition in diffusive systems for the barrier to be 
ineffective is £ 3> l n /F, where £ is the distance from the 
interface at which electrons and holes are still coherent. 
The coherence distance £ is determined by the length of 
the slab, energy of the electron, or magnetic field, de- 
pending on the physical case considered. When mea- 
suring the conductance of an NS junction, then for zero 
energy and zero magnetic field, the length of the nor- 
mal metal, d, is what limits the trajectories to lengths 
of order d 2 /l n (since a particle that reaches a distance 
d from the interface enters the reservoir, where phase 
coherence is lost). Therefore, the barrier is not effec- 
tive when r 3> l n /d. Since the transmission probability 
through the diffusive normal part is roughly T ss l n /d, 
this condition reduces to the known condition |j3"| for 
the barrier to be ineffective, T ^> T [in accordance with 
Eqs. @ and (§f)]. 

Following the same considerations one can obtain the 
condition for the barrier to be ineffective in various cases, 
noting the different mechanism impairing electron-hole 
coherence in each case. We consider, for example, the 
local DOS in a semi-infinite NS junction. The reduction 
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of the local DOS in the normal side of an NS interface 
is closely related to the averaged amplitude of an elec- 
tron near the interface to return to the same point as 
a hole through the pair amplitude (ipiipi) J2j|. At zero 
energy and assuming the normal metal is semi-infinite, 
there is no mechanism that limits the length of coherent 
trajectories. The electron hits the barrier as many times 
as needed, N ^> without losing electron-hole phase 
coherence, and according to Eq. (|J) , it is finally Andreev 
reflected. This results in a finite pair amplitude through- 
out the normal part, even in the presence of a barrier (at 
e = there is no reduction of the pair amplitude due to 
phase averaging), and in a zero DOS at zero energy. 

At finite energy e, the electron and hole moving in 
opposite directions in a trajectory of length L accumu- 
late a relative phase of Le/ ' (Tivf) which limits the length 
of coherent trajectories to order Tivf /e, and to distance 
£ = in = \JftD n /e from the interface (trajectories that 
traverse a distance longer than £„ result in phase differ- 
ence of order 2~k). The condition for having a large An- 
dreev reflection amplitude is therefore V ^> l n /£n [ [33[ . 
This condition assures that the total phase accumulated 
by the ingoing electron and outgoing hole is less than 2tt, 
and therefore the averaging in (V'xV't) results in finite pair 
amplitude and the local DOS is reduced. Similar consid- 
erations result in the width of the zero-bias anomaly in 
reflectionless tunneling being proportional to the Thou- 
less energy pjjf . 



VIII. CONCLUSIONS 

We have demonstrated the effect of reflectionless tun- 
neling in a ballistic NS system in which the multiple re- 
flections from the interface are due to the geometry. We 
considered a normal slab with superconductors attached 
to its sides, so the normal current flows in parallel to the 
NS interface. The barrier at the NS interface was taken 
to be smooth, so that normal reflection is specular, and 
with transmission probability r <C 1. We obtained a 
formula for the Andreev reflection amplitude from a tra- 
jectory that hits the barrier at the NS interface N times 
and showed that, when N 3> the barrier is trans- 

parent to pair current (though r < 1), leading to good 
proximity. 

We have shown that having a smooth rather than 
rough barrier at the interface is advantageous in giving 
rise to more pronounced and delicate features, which are 
not averaged over. This results in new measurable phe- 
nomena, such as the sharp peaks in the NS conductance 
as new channels open (in contrast to the usual step func- 
tion) and quasiperiodicity of the conductance as func- 
tion of magnetic field H. The smoothness of the barrier 
also enables one to conduct detailed manipulations such 
as extracting out a single channel from a normal metal 
(semiconductor) waveguide or extracting the current at 
a given position along the waveguide. 



By obtaining explicit formulas for the three-terminal 
conductances of the system as function of H and of the 
superconducting phase difference <& s , we have shown that 
both H and $ s impair the constructive interference lead- 
ing to the enhanced NS conductance, but in a qualitative 
different way. While as a function of H the enhanced 
NS conductance is limited to a small range of magnetic 
field and is channel specific, as a function of $ s the en- 
hanced NS conductance is generic, and is destroyed only 
near specific values of <£> s (odd multiples of it) for all the 
channels, leading to giant conductance oscillations. This 
difference is also reflected clearly in the shot noise behav- 
ior as function of both quantities. 

By demonstrating the possibility to obtain large An- 
dreev reflection in clean semiconductor-superconductor 
interfaces and the new possibilities such structures open, 
we hope to encourage experimental work in this regime. 

Our results were obtained using a semiclassical formal- 
ism, with which we reduced a two-dimensional nonsep- 
arable problem to an effective one-dimensional problem. 
We have demonstrated the usefulness of this formalism 
in a few situations, and hopefully it can be used in the 
future to solve other problems which are hard to tackle 
using the conventional techniques. 

We used this approach also for diffusive NS systems 
and demonstrated the connection between the effects of 
reflectionless tunneling in ballistic and diffusive NS junc- 
tions. We then considered the phenomena of reflection- 
less tunneling and the reduction in the density of states 
in diffusive NS junctions and showed that both can be 
obtained from a general criterion for the barrier, though 
large, to be transparent to pair current. 
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APPENDIX A: ANDREEV REFLECTION FROM 
AN N TRAJECTORY IN SNS JUNCTIONS 

In this appendix we describe the recursion formalism 
leading to Eq. ( |2p| ) and obtain a similar equation for tra- 
jectories that hit the NS interfaces an odd number of 
times. 

We choose the phase of the upper superconductor to 
be $ s /2 and the phase of the lower superconductor to be 
— $ s /2. The Andreev reflection amplitude of an electron 
hitting the upper (lower) boundary is 
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,+(-) 



(1) 



D +(-)t*./2 



(Al) 



and the Andreev reflection amplitude of an incoming hole 
is given by r e h(l) = — ?*£ e (l)- Throughout this appendix 
we consider trajectories that hit the NS interfaces any 
number of times, with the last hit occurring at the top 
interface. This is done for simplifying the calculation, 
and since the final result is even in <i> s , it does not depend 
on this assumption. However, the treatment for trajecto- 
ries that hit the interfaces an odd or an even number of 
times has to be done separately. For an odd trajectory, 
the recursion relation is 

r he (2N l)-r he (l)+ x _ _ 2) > (A2) 

where ^eeVhh)^) lii ^ e norma l reflection of an electron 
(hole) from the upper NS interface. For an even trajec- 
tory the recursion relation is 



r, (2N) r- (1) + ^(1)^(2^-1)^(1) 
rHe[ } hei)+ l-r- h (l)r he (2N-l) 



We define 



and 



r he (2N-l) = r^(l) 



P2N-1 

72JV-1 



2N 



12N 



(A3) 



(A4) 



(A5) 



Using the relations written after Eq. (||), we obtain 

02N-2 + 12N-2 



r he (2JV-l) = r+(l) 



02N-2 + 12N-2 



(A6) 



which we insert into Eq. (A3), and obtain the equation 



P2N 

72JV 



l + r^e- ltPs l + e 

r 2 (l + e i#s ) 1 + rV** 

.-<*« l + e -i*. 

e'*') 1 + rV* 8 



1 



: (1 



N 



P2N-2 
12N-2 



1 



(A7) 



where the last equation is obtained by explicitly finding 
that 



02 
72 



1 + e"**' 
1 + rV*' 



(A8) 



Following the same route for the odd case, we obtain 



P2N+1 

72W+1 



1 + r z e 



r 2 (l + e 

1 + rV* a 1 + e 
r 2 (l + e- l *0 l + r 2 e 



) l + r 2 e 



JV 



02N-1 
72JV-1 

1 



Diagonalizing the matrices in Eqs. (A7) and (A9) we ob- 
tain for the odd case 



R he (2N + !,$„)= \z + (l-z) coth 2 



y + N tanh 



2rcos($ s /2) 7(1^7) 
1 + r 2 cos $ s 



(A10) 



where 



y = tanh 1 



-i3>„/2 + ir 2 sm ($ s /2) 



(All) 



and for the even case Eq. (|20|), These equations are sim- 
ilar, only in the latter y = 0. 



APPENDIX B: VALIDITY OF THE 
CONDUCTANCE FORMULA FOR THE 
DIFFUSIVE SLAB 



In Sec. VII we obtain the linear conductance of a dif- 



fusive NIS junction using the approximation that the 
electron's motion in the normal slab is deterministic. 
However, for a normal slab with short range disorder 
Beenakker et al. p2| ] obtain a formula which differs in 
the T > T, T -> limit by a factor of 2 from Eq. (ph. 



In order to understand the factor of 2 difference be- 
tween the two cases we use the conductance formula for 
zero temperature of Beenakker et al. p2| 



4e 2 " T 2 



T V 



(Bl) 



where T m is the mth transmission eigenvalue of the nor- 
mal slab. 

The total transmission probability through the normal 
slab is given by T = Y,m =1 T m . In our approximation of 
the long-range scattering potential, each electron enter- 
ing the slab is predetermined, according to the position 
and direction at the entrance, to be either transmitted 
through the slab or reflected back to the normal reservoir 
(deterministic scattering). The transmission eigenvalues 
of the normal slab are therefore all zero and unity and in 
this case 



^ (2 

m=l v 



T 2 



T Y< 
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and 



4 e ^ " 4 e ^ 

771— 1 



(B3) 



increases, we have probability of order 1 to have a loop of 
length N 2 l, and therefore the length of the longest loop, 
as well as the length of the whole trajectory, is of order 
N 2 l. 



Since Gm — (2e 2 /h) Ylm=i ^mi we obtain the relation 
G NS = 2G N . However, in general, T 2 n / (2 - T m ) 2 < T m 
for all < T m < 1. In the case of short-range disorder the 
distribution of the transmission eigenvalues is such that 
£™ =1 T£/(2 - T m f = i£™ =1 T ro = \T H, and there- 
fore Gats = (2e 2 /h)T = Gn- This results in a factor of 
2 difference in the T — ► limit between our conductance 
formula ( p3f ) and the formula obtained by Beenakker et 
al. 



APPENDIX C: LENGTH OF A DIFFUSIVE N 
TRAJECTORY - RANDOM WALK THEORY 

We are interested in the question of how long a tra- 
jectory in the normal metal has to be in order to hit the 
interface N times. Looking at random walk in two di- 
mensions on a lattice which is rotated by 45° from the 
coordinate axes, it is easy to see that, since we are not 
interested at the exact point the trajectory hits the in- 
terface, there is a one-to-one correspondence between re- 
turning to the interface in two dimensions, and returning 
to the origin in the one-dimensional random walk. The 
question of return probabilities in one dimension is ad- 
dressed in Ref. Q, Chap. 3, using random walk path 
theory. This approach is elementary and very instruc- 
tive, and here we will just state its main results concern- 
ing our problem. Using path theory, Feller shows that in 
a one-dimensional random walk model, the probability of 
a first return to the origin after k steps is approximately, 
for large k, ff. — (2y / 7rfc 3 / 2 )~ 1 . Therefore, there is no 
average length for the first return (J2h=o^f k diverges). 
This peculiar result leads to a nonlinear dependence of 
the length of the trajectory on the number of times it hits 
the interface. (If there were an average return length a, 
then the average length of a trajectory that hits the in- 
terface N times would be aN). It is further shown that 
the probability to hit the interface iV times in a trajec- 
tory of length smaller than L is a function of L/N 2 l = w, 
where I is the mean free path and is given by 



P(w) 



-° 2 ' 2 ds 



(CI) 



This means that in order to hit the interface TV times 
a particle has to travel a length of order N 2 l. A typical 
trajectory of length N 2 l that hits the interface N times is 
not made of TV — 1 loops of similar length. The main con- 
tribution to the length of such a trajectory comes from 
one or two of its longest loops, whose lengths are of order 
N 2 l. This arises from the fact that: J2T=n 2 / fe ~ V-^j 
which means that if we have N returns to the origin, 
about one of them is going to be longer than N 2 l. As N 



[1 
[2. 
[3 
[4 
[5 
[6 

[T 

[8 
[9 
[10; 
[11 
[12 
[13 
[14 
[15 
[16 
[17; 

[is: 

[19 

[20 
[21 

[22 
[23 

[24 

[25 

[26 



A. Kastalsky, A. W. Kleinsasser, L. H. Greene, R. Bhat, 
and F. P. Milliken, Phys. Rev. Lett. 67, 3026 (1991). 

C. Nguyen, H. Kroemer, and E. L. Hu., Phys. Rev. Lett. 

69, 2847 (1992). 

B. J. van Wees, P. de Vries, P. Magnee, and T. M. Klap- 
wijk, Phys. Rev. Lett. 69, 510 (1992). 

J. Nitta, T. Akazaki, H. Takayanagi, and K. Arai, Phys. 
Rev. B 46, 14286 (1992). 

H. Takayanagi, T. Akazaki, and J. Nitta, Phys. Rev. B 
51, 1374 (1995). 

H. Takayanagi, T. Akazaki, and J. Nitta, Phys. Rev. Lett. 
75, 3533 (1995). 

A. Dimoulas, J. P. Heida, B. J. v. Wees, T. M. Klapwijk, 
W. v. d. Graaf, and G. Borghs, Phys. Rev. Lett. 74, 602 
(1995). 

A. F. Morpurgo, S. Holl, B. J. van Wees, T. M. Klapwijk, 

and G. Borghs, Phys. Rev. Lett. 78, 2636 (1997). 

A. F. Morpurgo, B. J. van Wees, T. M. Klapwijk, and G. 

Borghs, Phys. Rev. Lett. 79, 4010 (1997). 

J. P. Heida, B. J. van Wees, T. M. Klapwijk, and G. 

Borghs, Phys. Rev. B 60, 13135 (1999). 

P. A. M. Benistant, H. van Kempen, and P. Wyder, Phys. 

Rev. Lett. 51, 817 (1983). 

P. A. M. Benistant, A. P. van Gelder, H. van Kempen, 
and P. Wyder, Phys. Rev. B 32, 3351 (1985). 
J. H. Schon, C. Kloc, R. C. Haddon, and B. Batlogg, 
Science 288, 656 (2000). 

J. H. Schon, C. Kloc, and B. Batlogg, Nature 406, 702 
(2000). 

E. V. Sukhomkov and I. B. Levinson, Sov. Phys. JETP 

70, 782 (1990). 

A. Kadigrobov, A. Zagoskin, R. I. Shekhter, and M. Jon- 
son, Phys. Rev. B 52, R8662 (1995). 
N. R. Claughton, R. Raimondi, and C. J. Lambert, Phys. 
Rev. B 53, 9310 (1996). 

N. K. Allsopp, J. S. Canizares, R. Raimondi, and C. J. 

Lambert, J. Phys. Condens. Matter 8, L377 (1996). 

H. A. Blom, A. Kadigrobov, A. M. Zagoskin, R. I. 

Shekhter, and M. Jonson, Phys. Rev. B 57, 9995 (1998). 

N. M. Chtchelkatchev, JETP Lett. 73, 94 (2001). 

A. F. Morpurgo and F. Beltram, Phys. Rev. B 50, 1325 

(1994). 

C. Beenakker, Phys. Rev. B 46, 12841 (1992). 

M. Schechter, Ph.D. Thesis, The Weizmann Institute of 
Science, 2001. 

G. Wendin and V. S. Shumeiko, Superlattices Mi- 
crostruct. 20, 569 (1996). 

V. T. Petrashov, V. N. Antonov, P. Delsing, and R. Clae- 

son, Phys. Rev. Lett. 70, 347 (1993). 

V. T. Petrashov, V. N. Antonov, P. Delsing, and R. Clae- 



13 



son, Phys. Rev. Lett. 74, 5268 (1995). 
[27] M. P. Anantram and S. Datta, Phys. Rev. B 53, 16390 
(1996). 

[28] J. Torres, T. Martin, and G. B. Lesovik, Phys. Rev. B 

63, 134517 (2001). 
[29] A. L. Fauchere, G. B. Lesovik, and G. Blatter, Phys. Rev. 

B 58, 11177 (1998). 
[30] C. W. Beenakker and H. van Houten, Phys. Rev. B 43, 

12066 (1991). 



[31] We thank Nathan Argaman and Doron Cohen for clari- 
fying this point to us. 

[32] C. W. J. Beenakker, B. Rejaei, and J. A. Melsen, Phys. 
Rev. Lett. 72, 2470 (1994). 

[33] A. Volkov, Physica B 203, 267 (1994). 

[34] F. Zhou, B. Spivak, and A. Zyuzin, Phys. Rev. B 52, 
4467 (1995). 

[35] W. Feller, An Introduction to Probability Theory and Its 
Applications (Wiley, New York, 1966). 



14 



